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INTRODUCTION 


This report summarizes accomplishments of a Factory Mutual Research Corporation 
(FMRC) project on the Prediction of Fire Dynamics for the NIST grant period indicated. Work 
performed under a subcontract by Professor H.W. Emmons on Transient Ceiling Jet Heat Transfer 
is described as part of Task 1, Prediction of Fires in Buildings. This task has developed new 
practical expressions for heat transfer rates between the fire-induced ceiling-jet and the ceiling 
surface, so that the rate of increase in ceiling temperature can be predicted more reliably. 

‘The accomplishments of three tasks performed at FMRC are then presented in summaries 
of Tasks 2, 4, and 5. In Tasks 2 and 4, considerable progress has been made in the formulation 
of similarity relationships for gas temperature and soot layer thickness profiles in real-scale wall 
fires. We believe that these similarity relationships can be exploited to describe turbulent 
buoyant flames on vertical surfaces more accurately than ever. 

All of this work is aimed at the development of submodels or algorithms that can be used 
in NIST/BFRL comprehensive computer fire models. The FMRC tasks have previously led to 
the development of practical smoke-point measurement techniques for solid combustible _ 
materials, as described in a paper presented at the 4th International Symposium on Fire Safety 


Science, Ottawa, June 1994. 


TASK 1: Prediction of Fire in Buildings (H.W. Emmons, Gordon McKay Professor of 
Mechanical Engineering Emeritus, Abbott and James Lawrence Professor Engineering Emeritus) 
. Transient Ceiling Jet Heat Transfer 


A ceiling jet heat transfer-theory should be analyzed along with the dynamics$?At present a 
_ complete analyses do not agree as well with experiment as when the dynamics are assumed constant? In 
the following the dynamics are assumed to be known and are set at the mean values of velocity u, depth 6, 

~ and mass flow per unit width m. 
Under these conditions, the ceiling jet energy equations reduce to temperature equations 


describing the heat conduction and radiation from the ceiling jet to the ceiling and to ambient temperature 
objects below while the ceiling radiates to the ambient. 


The equations are: for the ceiling jet, 
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Although the open end of a corridor has no direct thermal effect on upstream heat transfer, it is a useful 
length scale for nondimensionalization. Thus introduce the variables 
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The corresponding local and average heat transfer rates are 
ceiling jet to ceiling 


"= h(To - T.8.- 8) 3 He 

Q= hy - TOL’@,- 6.) | , . iz 
ceiling jet to ambient below 
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ceiling, to ambient below 

Paneer je. a 

Q= halo - Tah 8, 17 
where the overbar indicates the arithmetic average over the effective length L” 

L’ =4L if ut>L 18 


x if x<L and ut>x 


Analysis of Results 


Case | An important case has a ceiling held at ambient temperature 6, = 0. The general solution 
reduces to 
6,= oe (Pst Pads U(ut-x) 19 
and the heat transfer rates ceiling jet to ceiling and ambient are: 
gq" =(h+ h,X(To - T.)8, 20 
Q=(h+hTo-TOL’O, = (h+hQL’AT 21 
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It is instructive to examine a specific case of ceiling, jet heat transfer in a corridor. The case 
chosen is as follows: Fire gases at 800 K issues at .5 kg/sec into a corridor 3m wide ( Mm = .1667kg/sec m) 
and length L = 100m. The ceiling is initially at 300 K but heats up with time. The ambient temperature 
below is at T, = 300 K and remains so. The properties of the ceiling, were taken as 1/2 ~ thick gypsum 
board for which p.d.cow = 1024 J/m*K . 

The mean ceiling jet properties were taken at T.y= 556 K with a density of p, = Pssu =.642kg/ a 
which implies an effective gravity of g? = 2Ap/p = 4.459m/sec”. Following Zukoski “, the ee jet 
volume flow per unit width Q=m/p, = .260m7/sec which implies a jet velocity of u = (g¢Q)" = 
1.050ny/sec and a jet depth of 8 = (Q’/g’)'? = 247m. Finally the connective heat transfer to *he ceiling, is 
given by Nu = .013 Re, Zukoski>. The thermal conductivity and viscosity of the ceiling jet was taken as 


that of air (mostly N2) at the mean film temperature of 425 K. Adding a linearized radiative heat transfer 
coefficient results in: 


h = 28.8w/m?K. which implies p, = .17]1 m! and Pw =.00293 nm! 
the remaining radiative heat transfer coefficients are 

h.=12.9 giving p, = .0763 nm 

h.=7.00 giving pw=.00055 m™ 


The solution for the ceiling jet temperature, 6,, equation 9 is presented as figure 1, and for the 
ceiling, 6, equation 10 as figure 2. The ceiling jet to ceiling temperature difference, @, - 8. equation 11 - 
figure 3, shows the local heat transfer rate from the ceiling jet to the ceiling, equation 12. As expected the 


€ 


heat transfer rate is large at the beginning, of the corridor and at early times and falls off approximately 
exponentially with time and location, Seanh Nera al ; payee 

For some purposes only the mean heat transfer rales are needed, equation 13, which is shown as 
the average temperature difference for length x and time ut in figure 4 for the specific case of a 10 meter 
long corridor. For this specific case it shows that the ceiling jet advances to the end of ts corridor in 
time t = 9.5sec after which the mean temperature difference in the corridor falls as the ceiling heats up 


Figure 5 shows that the mean heat transfer rate falls rapidly while the ceiling, jet first advances 
but falls slowly for a long time thereafter. . 


Since by plotting all of the temperature variations vs x on semilog, paper show only slightly bent 
lines, it would be expected that the frequently used log, mean of the end temperature differences would 
give a rather good simple approximation to the over all heat transfer. As the points x in figure 5 shows, 
the log mean approach gives heat transfers that are too high up to 42%. If experimental! data identical to 


the present theory were measured and were interpreted using, the log mean temperature difference, the 
maximum would be only about 20%. 


Drop from ceiling, 


ow 


The ceiling jet as it flows down a long, corridor has its buoyancy becomes so small that it cannot 
reinain at the ceiling against miscellaneous drafis, etc. The theory of this process is very complex and has 
not been adequately explored. To see the general effect, it was assumed that the ceiling jet could ses 
remain at the ceiling beyond 6,= .01, i.e. for the present case the ceiling jet temperature is within 5 K of 


the ambient. From figure 1 it is seen that the 8, = .01 point moves forward much more slowly than the 
normal jet velocity. This is shown in figure 6. 


Other heat transfer coefficients. 


The local heat transfer rate between the ceiling jet and the cciliny depends on the actual local 
temperature difference between them. .Since the actual ceiling jet temperature is given by a formula al 
least as complex as in this paper, it has been suggested that the ceiling jet heat transfer coefficient be 
based upon an adiabatic ceiling, jet temperature defined as the ceiling jet with no heat transfer to the 
ceiling. This “adiabatic” temperature is a simple exponential Cooper’ For the present case, the ceiling, jet 
actual temperature and the arbitrarily selected adiabatic femperature are compared in figure 6. The 
resultant heat transfer coefficients would be very different. For example, if the ceiling had a temperaturc 
distribution equal to the adiabatic ceiling Jet while the actual temperature of the ceiling jet were the lower 
line in figure 6, the adiabatic approach would be no heat transfer while in fact it is quite large. The use of 


the simpler adiabatic temperature can of course reproduce exactly the data from which it was obtained 
and is useful in fire calculations which are close to the conditions of the oriyvinal test. 


. 


1 Alpert R, Turbulent Ceiling Jet Induced by Large Scale Fires. Comb. Sci.Tech {1 p197 1975 
2 Emmons H.W. The Ceiling Jet in Fires NIST-GCR-90-102-1990 


3 Chobotov M.V. Gravity Currents with Heat Transfer Effects PhD theses CIT 1987 
4 Zukoski E. Personal Communication (Fr = 1) 
5 Zukoski E. Heat Tramsfer in Unwanted Fires ASME-JSME Thermal Eng. Conf. Honolulu 1987 


6. Cooper L. Woodlouse A The Buoyant Plume Driven Adiabatic Ceiling, Temperature Revisited. ASME 
Jour.Heat Trans. 108 p822 1986 


7 Abramowitz M. Stegun I Hdbk Math.-Func. NBS App Math. Series 55 1965. 
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TASK 2: Models for Turbulent Flame Chemistry and Radiation (J. de Ris) 
and 


TASK 4: Soot and CO Oxidation in Buoyant Turbulent Diffusion Flames (ie Orloff) 

The objective of Task 2 is to develop models for predicting 1) wall fire radiative feedback 
to the fuel surface, which controls wall fire burning rates, and 2) the combined soot and CO 
oxidation which controls the release of incomplete (toxic) products of combustion (soot, CO, etc.) 
and radiation from buoyant turbulent diffusion flames. Task 4 has the objective of providing 
experimental data for the development of these models. 

Considerable progress has been achieved during the past year towards the formulation of 
similarity relationships that may be very encouraging to fire modelers. Analysis of temperature 


’sl.2 radiant extinction measurements 


and soot layer thickness profiles together with Markstein 
show that the profiles remain geometrically similar for a fixed overall fuel to entrained air 
equivalence ratio corresponding to a fuel mass transfer rate, mn", which increases with the square 
root of height, Z. For additional detailed descriptions of the measurement and data-reduction 
techniques that were employed to reach this conclusion, see Appendix A-E. 

Various pyrolysis zone heights were simulated by supplying propylene to up to ten 132 
mm high and 320 mm wide water-cooled sintered-metal gas burners. The forward heat transfer 
zone was simulated by a 660 mm high water-cooled heat transfer plate mounted above the gas 
burners. The flow was maintained two-dimensional by 150 mm deep water-cooled side walls 
attached to the burner apparatus over its entire height. 

The thickness of the soot layer, 5,, was measured by inserting arrays of 5 mm glass rods 
into the flames perpendicular to the wall surface and rapidly withdrawing them after a two 
second exposure to the flames. The soot layer thickness was determined from the average length 
of soot deposit on ten rods. 

Figure 1 shows the length ratio, Z/6,, correlated against the inverse modified equivalence 
ratio p NIA */(th"-rn"). The modification of the equivalence ratio is required because the soot 
vanishes and the flames become blue for mass transfer rates less than mj=4 g/m’s. As a result 
of wall cooling and dilution by the products of combustion at low mass transfer rates, the data 
show that the cut-off, rh, is independent of Z. The turbulent motion typically transports the 
luminous soot (e.g. visible flames) out into the incoming air all the way to where the mean gas 


temperature drops to around 1000K. The correlation says that the soot layer thickness is 


proportional to Z for a given ratio of entrained air, p M@27Z Die to supplied fuel m"Z above its 
blue flame value mjZ. The correlation extends over a very wide range of flame equivalence 
ratios. 

Temperature profiles across the flame boundary layer were measured by a thermocouple 
rake consisting of 15 insulated Chromel-Alumel thermocouples inside 1.6 mm diameter Inconel 
sheaths spaced 12.6 mm on center and protruding 1 cm downward into the rising flow. The 
measured temperatures inside the flame were significantly depressed by radiation heat loss. On 
the other hand, the thermocouple temperatures outside the flame were significantly increased by 
radiant heat transfer from the flame. We corrected for both these effects with a simple heat 
transfer model to obtain the correlations shown in Figure 2. The vertical dashed line at y/d=1 
shows the boundary of the soot layer occurring at temperatures near T=1000K. Inside the soot 
layer, the presence of cold fuel is seen by the temperature drop near the wall. Outside the soot 
layer, the mixing of combustion products with the entrained air causes the corrected temperatures 
to asymptotically approach the ambient temperature as y/5 increases. These outer temperature 
profiles correlate reasonably well when plotted against (y-5)/Z. 

Markstein? reports measurements of the extinction of infrared radiation, £,, by soot (at 
wavelengths A, = 0.9 um and 1.0 um) across the flame boundary layer. These measurements 
immediately provide the integral of the soot volume fraction, f,6,, across the flame, f,5,=-A, 
In(1-€,)/7, which is replotted here in Figure 3 using coordinates similar to those in Figure 1. The 
factor of 7 in the above expression is recommended by Hottel and Sarofim?’ for soot. The 
correlation is not as good as for the soot standoff distance (Figure 2), apparently because the soot 
volume fraction, f,, depends weakly on the flow time. 

The faired curves in Figure 3 indicate that the fractional conversion of fuel carbon to soot, 
%X,» at a fixed mass transfer rate (1) initially increases proportional to the flow time up to a height 
Z=0.75 m for the C,H, flames, and then (2) becomes constant at greater heights. 

Figure 4 shows the measured transverse temperature profiles for four C,H, flames at two 


heights and two equivalence ratios, rn/Z” near 11 and 26 g/m?” 


s respectively. Notice the 
dependence of profile-shape on equivalence ratio, but almost perfect similarity at fixed 
equivalence ratios. Theory suggests that the heat release by combustion per unit wall height 
increases with Z’” at fixed equivalence ratio. The present data suggest that the soot and gas 


radiation initially increase almost linearly with height for small Z where radiation is unimportant, 


but convective heat loss is important. At somewhat greater heights the radiation increases more 
nearly with Z'? due to radiant heat loss. Ultimately the radiation will level off at very large Z 
after the flames become optically thick. 

The similarity of combustion observed in this study greatly eases the task of analyzing 
experimental data and will make it much easier to tailor future detailed semi-empirical turbulent 
combustion models, so that their predictions exactly agree with experiment (at least for situations 
having m" ~ Zry. One rarely finds such true similarity (i.e. reducible to one-dimensional 
behavior) in combustion, especially turbulent buoyant combustion. 

We believe that these similarity relationships can be exploited to describe turbulent 
buoyant wall fires, with the radiative and completeness of combustion properties of the fuel 
characterized by the smoke point. A series of radiant and total heat feedback measurements in 
turbulent wall fires are planned in which the use of binary C,H,/C3H, fuel mixtures will provide 
comparisons of these measurements with predictions based on the smoke point of the fuel 
mixture. Markstein’s~ measurements on C,H, wall fires have shown that the overall heat transfer 
back to the surface is only about 70% of the outward heat transfer. This important effect of 
radiant blockage due to cold sooty gases near the wall will be studied in the planned series of 


C,H,/C3Hg fires. 


References 

1. Markstein, G.H. and de Ris, J.: Twenty-Fourth Symposium (International) on Combustion, 
1747, The Combustion Institute, 1992. 

2. Markstein, G.H.: Personal Communication, 1994. 

3. Hottel, H.C., and Sarofim, A.F. (1962), Radiative transfer, McGraw-Hill, New York, p. 199. 
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Task 5: Methods for Developing Probability Estimates in Fire Risk Analyses (D.M. Karydas) 
The Fire Risk Analysis approach establishes a performance based methodology where the 
performance parameter is quantified risk associated with fire hazard. The general framework of 
the approach includes the deterministic models of fire spread and smoke movement, the 
probabilistic models of the safety and control systems, and the event tree models with fire 
initiating event and engineered safety system functions. The associated concepts of accident 
prevention, consequence mitigation and hazard elimination are established. Inputs to the event 
tree models are provided by the deterministic models of fire spread and growth, smoke movement 
and deposition, and the detection/protection actuation. Additionally, the frequency of fire 
initiating and the success/failure probabilities of the fire safety and control systems (such as 
detection systems, protection systems, etc.) are seoruks by the probabilistic models. Adjustment 
to site specific fire risk assessment that deviate from the generic models are introduced at the 
level of accident initiation or consequence mitigation by using Bayesian inference theory. Based 
on this fire risk framework, a computer program is under development that can demonstrate the 
functional relationships of the elements of the methodology outlined and provide on-line 


explanation of the basic fire risk assessment terminology. 
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APPENDIX A 
SCALING OF SOOT AND TEMPERATURES IN WALL FIRES 
L. Orloff 


Introduction 
The present study is directed towards the measurement and analyses of optical path length, 
temperature profiles, and soot layer thickness in turbulent wall fires. This work parallels 


Markstein’s! 


investigation of the radiant heat transfer blockage by cold gas and soot near the fuel 
surface. The results of these two studies have yielded new insights relating to the wall fire 
boundary layer, air/fuel mixing processes, radiative heat feedback to the fuel surface, convective 
heat release, and soot distribution and re-radiation in turbulent wall fires. 

In earlier work, upward fire spread and steady burning along vertical PMMA slabs was 
investigated”, and gaseous-fuel diffusion flames on flat surfaces of variable orientation, 
including the vertical, have been studied.* The role of radiative energy transfer was an essential 
element in understanding wall burning in these studies. Theoretical models of wall burning have 


79 However, all these models lack a 


also been developed at this laboratory>* and elsewhere. 
predictive input for flame radiation, which controls the burning process. We believe that 
Markstein’s and our results will provide some of the experimental data needed for the 
development of submodels (i.e., formulae or correlations) of wall heat transfer components for 


inclusion in upward spread and fire growth models. 


Description of the Experiment 
The pyrolysis region of a wall fire was simulated by supplying propylene to a vertical 
sintered porous metal surface. The burner consisted of ten 132-mm high and 380-mm wide 
panels topped by a 660-mm high solid metal heat transfer plate. To approach two-dimensional 
flow conditions as closely as possible, water cooled side walls of 150-mm depth were attached 
to the burner over its entire height. Further details of this apparatus are given in Markstein’. 
The fuel supply to the burner was metered by a mass flow controller and measured with 


an inclined manometer providing the differential pressure across the orifice of a laminar flow 
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element. The burner was situated in a water-cooled enclosure in a laboratory supplied with a 
quiescent flow of air. Product gases were vented through a passive hood directly above the 
burner. Fuel mass transfer rates, mh", up to 60 g/m’s could be obtained with three burner units 
operating (total height 366 mm). Activating additional bummer units reduced the maximum rm" 
due to a fixed maximum fuel supply rate. 

The soot layer thickness was measured by inserting arrays of 5 mm diameter glass rods 
into the flame perpendicular to the wall surface and rapidly withdrawing them after a 2 second 
exposure to the flames. Magrides and Dobbins!° have shown that the amount of soot deposited 
is proportional to the local soot volume fraction in the gas phase. Two arrays, each containing 
five rods were inserted for each measurement. The soot layer thickness was determined from the 
average length of the soot deposit on the ten rods. The 2 second exposure was chosen because 
shorter exposures sometimes failed to provide a well-defined soot trail, and longer exposures 
tended to leave thick soot deposits with poorly defined tails. The soot layer thickness is 
measured from the end of the rod placed against the wall to the location where the soot deposit 
is visually judged to have decreased to 50% of the maximum deposit on that rod. Often this 
position is manifested by a distinct demarcation between dark and lighter soot deposits. 

The rods were laterally spaced 41 mm apart and the array was centered on the midpoint 
of the 381 mm wide burner surface. Typically the standard deviation of ten soot lengths was 
10% of their average length. Soot layer thickness measurements were obtained for CH, supply 


rates, tn", 4 g/m’s - 60 g/m’s at heights, Z, 345 mm-1138 mm above the leading edge. 


Experimental Results 

Fig 1 shows the variation of soot layer thickness,6,, with mass transfer rate at four 
heights. 

The thickness,6,,approaches zero at all flame heights as the mass transfer rate approaches 
4 g/m*s. This critical mass transfer rate corresponds to a B-number of around 0.5, at which the 
flames turn blue due to a combination of wall cooling and dilution of the supplied fuel by 
products of combustion diffusing back to the wall. For mass transfer rates less than 4 g/m’s a 
brown liquid rather than soot is deposited onto the glass rods. 

Figure 2 shows the length ratio, Z/6,, correlated against the inverse modified equivalence 


ratio p MagZ yu */("-th,"). The modification of the equivalence ratio is required because the soot 


13 


vanishes and the flames become blue for mass transfer rates less than rh," = 4 g/m’s. As a result 
of wall cooling and dilution by the products of combustion at low mass transfer rates, the data 
show that the cut-off, rh, is independent of z. The turbulent motion typically transports the 
luminous soot (e.g. visible flames) out into the incoming air all the way to where the mean gas 
temperature drops to around 1000 K. The temperature of the luminous soot at this location 
appears to be much hotter. The correlation says that the soot layer thickness is proportional to 
Z for a given ratio of entrained air, p Rez) * to supplied fuel r"Z above its blue flame value 
m,"Z. The correlation extends over a very wide range of flame equivalence ratios. The turbulent 
gas motion causes the soot to diffuse outward well into the oxygen-rich region where it 
encounters the entrained ambient air being drawn towards the wall. 

Temperature profiles across the flame boundary layer were measured by a thermocouple 
rake consisting of 15 insulated Chromel-Alumel thermocouples inside 1.6 mm diameter Inconel 
sheaths spaced 12.6 mm on center and protruding 1 cm downward into the rising flow. The rake 
holder contained thermocouple connections to a fast scanning data acquisition unit. The rake was 
angled 56° from the normal to the burner surface for most measurements. This orientation 
provided 7.1 mm thermocouple spacing normal to the wall. Figures 3a and 3b show the 
measured thermocouple temperatures at two different heights over a wide range of mass transfer 
rates. Measured temperatures were averaged over 60 scans taking 95 seconds. The temperatures 
were not influenced by soot deposits over this measurement period. 

The measured temperatures inside the flame were significantly depressed by radiation heat 
loss. On the other hand, the thermocouple temperatures outside the flame were significantly 
increased by radiant heat transfer from the flame. Using knowledge of flame radiation (reported 
previously) we corrected for both these effects with a simple heat transfer model to obtain the 
correlations shown in Figures 4a and 4b. The model is described in Appendix D. 

Figures 4a and 4b show that the temperatures inside the soot zone are well correlated 
when the distance from the wall is normalized by the soot stand-off distance. The vertical dashed 
line at y/5=1 shows the boundary of the soot layer occurring at temperatures near T=1000K. 
Inside the soot layer, the presence of cold fuel is seen by the temperature drop near the wall. 
Outside the soot layer, the mixing of combustion products with the entrained air causes the 
corrected temperatures to asymptotically approach the ambient temperature as y/6 increases. 


These outer temperature profiles correlate reasonably well when plotted against (y-5)/Z. 
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Markstein! reports measurements of the extinction, €,,by soot of infrared radiation (at 
wavelengths A, = 0.9 um and 1.0 ym) across the flame boundary layer. His extinction results 
scaled by th"Z'? are shown in Appendix E. Curve fits to these measurements immediately 
provide the integral of the soot volume fraction, f,5., across the flame, 


Ey EA) 


£,5,=-A, In(1-€,)/7 


which is replotted here in Figure 5 using coordinates similar to those in Figure 2. The factor of 
7 in the above expression is recommended by Hottel and Sarofim!! for soot. The correlation is 
not as good as for the soot standoff distance (Figure 2), apparently because the soot volume 
fraction, f,, depends weakly on the flow time. 


The faired curves in Figure 5 indicate that the fractional conversion of fuel carbon to soot, 


X,» at a fixed mass transfer rate (1) initially increases proportional to the flow time, V 2279 , up 


to a height of 0.75 m for the C,H, flames, and then (2) becomes constant at greater heights. 
Taking the ratio of the ordinates of Figures 2 and 5 one evaluates the soot volume 
fraction, f,, shown in Figure 6 for each of Markstein’s C,H, wall flames. These volume fractions 
have considerable scatter, but are reasonably independent of flow conditions. The soot volume 
fractions around f, = 10° provide radiant heat fluxes in general agreement with those measured 


in C,H, flames. 


Combustion Similarity for m" (z) ~ Zo3 

Both the equations for turbulent buoyant combustion of boundary-later diffusion flames 
and measured transverse profiles of velocity’, temperature’® and species concentrations’, all show 
almost perfect similarity at different heights, z, in the pyrolysis zone when the supplied fuel mass 
transfer rate, mm" (z) increases with the square root of z. For such similarity flows: (1) all macro 
velocities are proportional to z”; (2) all transverse positions are proportional to z; while (3) all 
macro temperatures, densities and species concentrations are invariant at equivalent homogenous 
positions. It follows that the overall chemical heat release, q,, per unit area across the flames 


is proportional to z” and can be expressed by 
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Ga (z) = 21/2 £() (1) 


where f() is some function of the overall fuel to entrained air equivalence ratio of the flame 
vA Z 

) =f Tia, / (So vjieoz (2) 
@) Oo 


which is determined by the value of m" (z)/z” for a particular similarity flame. Figure 7 shows 
the measured transverse temperature profiles for two sets of C,H, flames, at two different 
equivalence ratios. Notice the dependence of profile-shape on equivalence ratio, but almost 
perfect similarity of profiles for fixed equivalence ratio determined by m" (z)/ z*. Soot and 
gaseous flame radiation deviate slightly from similarity, increasing at first almost linearly with 
height for small z, then increasing closer to z” as a result of radiant heat loss from the flames 
, and finally leveling off at very large z after the flames become optically thick. The convective 
heat loss to the wall has an opposite trend, being at first proportionately moderately large for 
small z, and then soon becoming negligible for large z as a result of mass transfer blowing. The 
combined radiative and convective flame heat transfer for the pyrolysis zone beyond z = 100 mm 
should be roughly similar (e.g., roughly proportional to z”) for hazardous-scale fires. 

This similarity of combustion greatly eases the task of analyzing experimental data and 
will make it much easier to tailor future detailed semi-empirical turbulent combustion models, 
so that their predictions exactly agree with experiment (at least for situations having 
tn" ~ (z)/z”). This is, so far, the only buoyant turbulent combustion situation which has been 
shown to be truly similar (e.g. reducible to a one-dimensional problem). 

More recently we have been exploring implications of theoretical arguments which 
suggest that the temperature and species concentration profiles across the boundary layer should 


be correlated by y/Z at fixed values of the flame equivalence ratio m"/Z’”. 


Conclusions 
1. The soot layer thickness, 6, is proportional to height, Z, for a given ratio of entrained air 


to supplied air above its blue flame value, m,"Z. The correlation extends over a very wide range 


of flame equivalence ratios. 
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2. Longitudinal flame temperature profiles inside the soot layer obtained for a range of fuel 


mass transfer rates are well correlated when the distance from the wall is normalized by the soot 


layer thickness. 


as The temperature at the boundary of the soot layer was 1000-1100 K for the C,H, flames 
studied. 


4. The fractional conversion of fuel carbon to soot, y;, at a fixed mass transfer rate (1) 


initially increases proportional to the flow time, V22/9 , up to a height of 0.75 m for the CH, 


flames, and then (2) becomes constant at greater heights. We hypothesize that soot formation is 
proportional to available flow time for the C,H, flames below 0.75 m. Above this height the 


available fuel carbon is increasingly oxidized. 


5. Soot volume fractions, f,, around 10° were obtained from the present 5, measurements 


and Markstein’s’ extinction data. Calculated radiant heat fluxes based on this finding are in 


general agreement with Markstein’s radiant heat transfer measurements. 


6. Measured transverse temperature profiles show nearly perfect similarity of profiles for 
fixed equivalence ratio determined by m" (z)/ z*. This discovery of a buoyant turbulent 
combustion similarity condition greatly simplifies the development of semi-empirical turbulent 


combustion models. 
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Figure 1. Variation of soot layer thickness with fuel mass flux at four heights above the leading 
edge. 
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Figure 2. Height to soot layer thickness ratio correlated against inverse modified equivalence 
ratio. 


20 


1000 


Moasurod Thermocouple Temperatura, % 


8 


8 


3 


8 


8 


3h Wall Fire 
TC tree atZ=365mm 


490 so 69 7a 


80 


Distance from Wall, Y (mm) 


1000 


re 
Moagured Thermocouple Tomperature, -¢ 


g 


3 


8 


3 


8 


8 


8 


100 6110) 6120 


21 


192 KW 293 g/sm* 
©" 163 KW 24.8 gsm 
"144 KW 220 g/sm~ 
"125 KW 19 gis m* 
—~ 107 KW 163.g/m~*s 
~~ SE KW 13.1 gim“s 
->"66.2kKW 10.1 gim~s 
"45.9 kW 7.0 gim*s 


- 26 KW 3.98 g/sm* 


—-16.4kW 25 gsm” 


122 kW 1.86 gins 


(3a) 


Cah, Wall Fire, 
TC tree at Z = 1022 mm 


Distance from Wall, Y (mm) 


375 KW 171 gism* 


-©- 330 KW 1531 g/sm~ 
-£-280 KW 128 g/sm~ 
—--273 KW 102 gism* 
——176 kW 8.0 gms 
~=-116 KW S3.gim7s 


"94.9 KW 4334 gim*s 


(3b) 


Figure 3. Measured thermocouple temperatures at heights (a) Z = 373 mm and (b) Z= 1022 mm 
for a range of fuel mass transfer rates. 
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Figure 4. Thermocouple temperatures of Figs. 3a and 3b respectively corrected for radiation heat 
loss by thermocouples inside the flames and radiative heat flux to thermocouples 
outside the flames. Vertical dashed line shows boundary of soot layer. 
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Figure 5. Ratio of height to integral of soot volume fraction across the flame correlated against 
inverse modified equivalence ratio. 
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Figure 6. Soot volume fractions in C,H, wall fires. 
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APPENDIX C 


CORRECTED THERMOCOUPLE MEASUREMENTS 
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APPENDIX D 


| EFFECTS OF FLAME RADIATION AND SURFACE RERADIATION ON 
THERMOCOUPLE GAS TEMPERATURE MEASUREMENTS IN AND NEAR FLAMES 


Thermocouples are often used to measure local gas temperatures both inside and near to 
flames. However, the radiative heat transfer from the flames increases the thermocouple 
temperature, while the radiative heat loss from the hot thermocouple surface decreases its 
temperature. In steady state, the convective heat transfer to the thermocouple from the gas 
flowing over the thermocouple must equal the radiative heat loss from its surface minus the rate 


of heat gain from nearby flames, or | 
DA GeagT) sou) = To)eke cle 7.) 


per unit thermocouple surface area. Here, T, T,, and T;, respectively, are the thermocouple, gas, 
and effective flame radiation temperatures; h is the convective heat transfer coefficient; o is the 
Stefan-Boltzmann constant 5.7 W/cm? (1000K)*; and €, 1s the average flame emissivity, as 
viewed by the thermocouple. We implicitly assume that the thermocouple is coated by a thin 
layer of black soot, causing it to have unit surface absorptivity and emissivity. Solving for T,, 


one has 


Te ee TOS O/ PD (4s eTe) 4G) EG elas le 


Convective Heat Transfer 
Our thermocouples point downward into the rising hot flow. Eckert and Drake (1972) 
provides the convective heat transfer formulas: 
Nu = hd/k = [0.43 + 0.50 Re®*]Pr°*® (cylinders in cross flow) 
Nu = hd/k = 2.00 + 0.236 Re®®°pr° (spheres) 


Both these formulas yield a Nusselt number around 4 for a Reynolds number 
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R, = ud/v = ¥2gZ d/v= 45 


Evaluated for a d = 0.15 cm diameter thermocouple in a one meter high (Z = 1m) flame at 1300 
K. The Reynolds number is quite uniform across the flame because the temperature dependence 
of the kinematic viscosity, v, causes it to decrease almost proportionally to the velocity p1, as the 
velocity decreases from its maximum value in the flame to its inward directed entrainment value 
in the cold gas surrounding the flame. Thus, modeling the thermocouple by either a cylinder or 


sphere, one has 


tee Niles Caseae eet) Ts) -2/ de= 001 (a) 1.) °°" W/cm K 


which takes into account the temperature dependence of the gas thermal conductivity, k. 


Flame Radiances 

Consider a planar flame of uniform thickness, 6, adjacent to a vertical wall. It is assumed 
that the effective flame radiation temperature, T;, is independent of the distance, y, from the wall, 
while the probability of finding a flame at position y decreases with y, so that the time-averaged 
flame absorption-emission coefficient a (y) decreases with y and becomes zero for y > é. 
Defining the local optical depth normal to the wall 


n= [7 @ dy 


oO 


The flame radiances along rays at an angle 6 w.r.t. the normal become 


N, (1,0) = 2 (Tears) (lee 16?) outward 


Nene Oia Leet) fee taal acco inward 
T 
where 7),, is the maximum value of 7 corresponding to y = 5. The outward flame radiance, Ne 
continues constant along a ray extending beyond the flame thickness y > 6. It only depends on 


the angle @. The inward flame radiance occurs only inside the flame, y < 6. 


Radiation to a Thermocouple Inside the Flame 
In general, the radiant heat transfer, q", impinging on a surface, is- 


q" =JN(o) coso dQ 


where Q is the solid angle and @ is the angle of a ray to the outward normal. In the case of a 
planar flame between two parallel surfaces, the integral over @ can be approximated by 
evaluating N(@ cos 6 at its "mean beam" value (@ = 7/3) and multiplying it by the solid angle 
Q = 2n facing the flame, to obtain 


yay al ats} 4) 
q! = N()cos@- O = © (rt - Tr‘) 428) 2n=0 (Te - To) tee 


a 


where 71) is the flame optical depth normal to the surface. The approximation is slightly improved 
by setting seco = 1.8 instead of 2 in the exponent. 

To evaluate the heat transfer to a spherical thermocouple of diameter d located inside the 
flame, we calculate the total inward and outward radiation impinging on a thin circular disk of 


area 1d7/4, parallel to the wall, to obtain the average heat transfer, q", per unit thermocouple area 


2 fe ek, re q- ae Ta /4 4 4 = 8 SA Oka yD 
qu = —— = — 6 (T; - Tz) [(1-e%'* 9) + (1-e Ly 
md? 1d? (72 ) ais 


This radiation is maximum at the center of the flame y = 1),, /2 and is minimum at both the inner 


edges y = 0, 7),;: 


a. = $ (Tr & 7 (1-e 70-2") 


CRE ce alia ty bee 


which become equal for optically thin flames n,, > O. 
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Radiation to a Thermocouple Outside the Flame 


For y > 6 the radiances emitted by the flame remain constant along each ray. However, 
the lateral side walls limit the extent of the flame and reduce the total solid angle of the flame 
as viewed by the thermocouple. Consider a uniform flame of thickness 6 in the form of a 
circular slab of radius R viewed by a thermocouple disk of area md7/4 located a distance y along 
the axis from the flame. Assume 6/y and d/y are very small. 


The total solid angle of the flame, Q, as viewed from y is 
8, . im 
O,=2n /[* sinOc® = 2m (1-cose, ). 
Oo 
The angle of the "mean beam" 6, subtends half this solid angle, 
Oye eee C050, ayer (1 -cosg.) 


so that 
ges0 e-leeecOsun)) 2 (iat y/yy-tn* )*/2 


Evaluating N(o) coso along the "mean beam" the heat transfer to the thermocouple (i.e., disk) 


per unit thermocouple area becomes 


2 
: 4 
ee note N(6,) cos8, Q, 
= = (TT) (1-e° 2a8e°®2) cos, -2n (1-cos6,) 


uO 4 4 0.9n,sec8, R?2 
ae - Te) (1- ———— 
4 (Tr ) ( S ) (vie +R?) 


where 


SECO y eee ey ly Vice R2) 
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After letting R go to infinity, this heat transfer smoothly matches the inside flame value at the 


flame edge y = 6. 


Flame Radiation Inside and Outside the Flame 


The above results can be combined into a single expression 


Tre R? 
oid es = (T? ws Te) a 


=| f(y) 


f(y) u (tee ate ces + h(6-y) ae Cee ry 
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APPENDIX E 


MARKSTEIN’S SOOT EXTINCTION MEASUREMENTS 


Markstein’s! 


measurements of the extinction, ¢, by soot of infrared radiation (at 
wavelengths 4, = 0.9 pm and 1.0 um) across the flame boundary layer shows that soot optical 
depth scaled by m"Z'”” is curve-fitted by [A/(50 + m")] {(th"+1h,")? - ("+ 4)*]!? where 
mn," = 4 e/m-s. and A= 0.6, 1.35 and 1.7 at Z = 345 m, ‘742, and 1-7 at Z = 0.345 m, 0.742 


m, and 1.138 m, respectively. 
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